Replicas of the Fano resonances induced by phonons in a subgap Andreev tunneling 
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We study influence of the phonon modes on a subgap spectrum and Andreev conductance for 
the double quantum dot verticaUy coupled between a metaUic and superconducting lead. For the 
monochromatic phonon reservoir we obtain replicas of the interferometric Fano-type structures 
appearing simultaneously in the particle and hole channels. We furthermore confront the induced 
on-dot pairing with the electron correlations and investigate how the phonon modes afi^ect the zero- 
bias signature of the Kondo effect in Andreev conductance. 
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I. INTRODUCTION 

Electron transport through nano-size transistors con- 
taining the quantum dots, molecules and/or nanowires 
is determined by the available energy levels (tunable 
by external gate voltage) and strongly depends on the 
Coulomb interactions Discretization of the energy 
levels is responsible for oscillations of the differential con- 
ductance upon varying the gate voltage, whereas the cor- 
relation effects lead to the Coulomb blockade and can 
induce (at low temperatures) the Kondo resonance en- 
hancing the zero-bias conductance to a unitary value — 
Besides promising perspectives for the applications 
in modern electronics/spintronics the nanoscopic struc- 
tures represent also valuable testing grounds for probing 
the many-body effects. Magnetic, superconducting or 
other types of orderings absorbed from the external leads 
can be confronted with the on-dot electron correlations 
in a fully controllable manner. 

In this regard, especially interesting are the hetero- 
junctions where the quantum dots (QDs) are in con- 
tact with the superconducting (S) electrodes. Nonequi- 
librium charge transport can occur there either via the 
usual single particle tunneling (upon breaking the elec- 
tron pairs) or by activating the anomalous (Andreev or 
Josephson) channels. The resulting currents are sensitive 
to a competition between the induced on-dot pairing and 
the Coulomb repulsion. In such context there have been 
experimentally explored the signatures of 7r-junction 
Josephson effect [5|, superconducting quantum interfer- 
ence M, quantum entanglement by splitting the Copper 
pairs [T|, multiple Andreev scattering 3: interplay 
of the on-dot pairing with the Kondo effect [ol-flTI . These 
and similar related activities have been discussed theo- 
retically by various groups [T2| - [T6| . 

Since in practical realizations the nanoscopic objects 
are never entirely separated from an environment (e.g. 
a given substrate or external photon/phonon quanta) 
therefore transport properties are also affected by the 
interference effects. A convenient prototype for studying 
such phenomena is the tunneling setup shown in figure 
[U where the central quantum dot is coupled to the side- 
attached quantum dot and eventually to other degrees of 
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FIG. 1: (color online) Schematic view of the double quantum 
dot coupled in T-shape configuration between the metallic (N) 
and superconducting (S) electrodes where an external phonon 
bath affects the side-attached quantum dot (QD2). 



freedom. Transport through this T-shape double quan- 
tum dot (DQD) occurs predominantly via the central 
quantum dot (QDi), whereas electron leakage to/from 
the side-coupled quantum dot (QD2) brings in the inter- 
ference effects. In the case of both normal (N) electrodes 
and assuming a weak interdot coupling it has been ar- 
gued [13, [13 that the differential conductance should re- 
veal the asymmetric Fano-type lineshapes. This fact has 
been indeed observed experimentally [19| . 

Recently we have explored the interferometric patterns 
for the DQD case placed between the metallic and super- 
conducting electrodes [20|. In such heterostructures the 
interferometric lineshapes appear simultaneously at neg- 
ative and positive energies because of the mixed particle 
and hole degrees in the effective quantum dot spectrum 
[2l| . These effects manifest themselves in the subgap 
Andreev conductance [l^l- Stability of interferometric 
Fano structures on a dephasing by external fermionic 
bath have been analyzed in Refs [22,, ^23,] . Here we ex- 
tend such study addressing the role of bosonic bath in 
the setup displayed in Fig. [TJ We argue that monochro- 
matic phonon bath induces a number of the Fano-type 
replicas depending on the adiabadicity ratio A/wq. 
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In the next section we introduce the microscopic model 
and briefly specify characteristic energy scales. We also 
discuss the formal aspects concerning the adopted ap- 
proximations. In the section III we analyze spectroscopic 
fingerprints of the phonon modes for the case of uncor- 
rected quantum dots. Finally, in section IV, we address 
the correlation effects (Coulomb blockade and Kondo 
physics) due to the on-dot repulsion between the opposite 
spin electrons. Appendix A provides phenomenological 
arguments for the Fano-type interferometric patters of 
the double quantum dot structures. 



II. MICROSCOPIC MODEL 

The double quantum dot heterojunction shown in Fig. 
[1] can be described by the Anderson- type Hamiltonian 



H 



DQD 



(1) 



where Hieads = Hjq + Hs denote the normal N and su- 
perconducting S charge reservoirs, Hjjqjj refers to both 
quantum dots (together with the phonon bath) and the 
last term Ht describes the hybridization to external 
leads. We treat the conducting lead as a Fermi gas 
Hn = X^k cr ^kAfcj^o-JV^ko-AT and we assume that isotropic 
superconductor is described by the BCS Hamiltonian 

Hs = Ek,a fkscLsCk<TS - A Ek(clct5clk;s + c-kisckts) • 
The operators oj^lp correspond to annihilation (creation) 
of the itinerant electrons with spin a i and energies 
Ck^a =£k^— MyS ai"^ measured with respect to the chemical 
potentials 

The double quantum dot along with the phonon bath 
is described by following local part 

HdQD = X! ^^<^^ad^a + t ^ (d\j2a+^-C^ (2) 

-I- ^ [/i rfl-t-^^it '^ii'^ii + ^oa^a + X! 4cr'^2(T(a^ + a), 

where we use the standard notation for the annihilation 
(creation) operators (^'^ for electrons at each quantum 
dot QDi=i,2. Their energy levels are denoted by ei and 
Vi refer to the on-dot Coulomb potentials. Since we are 
interested in the Fano-type interference we focus on the 
electron transport only via the central quantum dot 



I3=N,S k,CT 



(3) 



This situation can be extended to more general cases 
when electron tunneling directly involves both the quan- 
tum dots. For clarity reasons we postpone such analysis 
for the future studies. 



Gi{Ti,T2) = —iTr{^i{Ti)'^l{T2) defined in a representa- 
tion of the Nambu spinors ^| ee {dl^,di^), = (*|)'''. 
In the equilibrium conditions (for /i^v = Ms) such ma- 
trix Green's function depends only on time difference 
Ti— T2. The corresponding Fourier transform can be then 
expressed by the Dyson equation 



G-i(..)=fif-i(..)-SO(c.)-Sf(a;), 



(4) 



where the bare propagators gi{uj) of uncorrelated quan- 
tum dots are given by 



OJ - £i 






uj + 



(5) 



First part of the selfenergy (w) comes from a com- 
bined effect of the interdot coupling, hybridization of the 
central QDi with the external leads ([3]) and the phonon 
bath contribution acting on QD2. The other term S^(cj) 
appearing in Q accounts for the many-body effects orig- 
inating from the on-dot Coulomb repulsion Ui. 

Let us begin by first specifying the selfenergy S^C'^) 
for the uncorrelated central quantum dot. The usual di- 
agrammatic approach yields 



k,/3 



(6) 



where g^(k, w) denote the matrix Green's functions of 
the leads. In particular, we have for the normal lead 



g^(k,c.)= ( 



and for the superconducting electrode 



(7) 



gs(k,w) = 



"k I "k -"ktlk I Mkfk \ 

_ I w-Bk i^+Ek w-Bk w+Bk (g) 
-Mk-"k I "k"k "k I ^k / 



with quasiparticle energy E'k Vsks + A^ ^^"^ BCS 



coefficients 



'k' "^k 



wide band limit approximation we assume the constant 
hybridization couplings Tp — 27rEk l^k^P <5(w— Ck/s) and 
treat Fjv as a convenient unit for energies. We then for- 
mally have 



N 



k 

E 

k 

where 03 
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(9) 



(10) 



Outline of the formalism 



Energy spectrum and transport properties of the sys- 
tem ([T]) can be inferred from the matrix Green's function 



7(0;) 



for |a;| > A, 
for Iwl < A. 



(11) 



Deep in a subgap regime (i.e. for \uj\ ^ A) only the off- 
diagonal terms of (|10p survive, approaching the static 



3 



value — r5/2. From the physical point of view a magni- 
tude |— r5/2| = Adi can be interpreted as on-dot pairing 
gap induced in the QDi. Such situation has been stud- 
ied in the literature by a number of authors [13] apply- 
ing various methods to account for the correlation effects 



B. Influence of the phonon modes 

Generally speaking, the phonon modes have both the 
quantitative and qualitative influence on the electron 
transport through nanodevices [25|. In particular they 
can be responsible for such effects as: appearance of the 
multiple side-peaks, polaronic shift in the energy levels, 
lowering of the on-dot potential Ui (even to the nega- 
tive values promoting the pair hopping [1^), suppression 
of the hybridization couplings Tp and often serve as a 
source of the decoherence. These and related subjects 
have been so far studied by many groups, mainly con- 
sideriiig the single quantum dots coupled to the normal 
leads [27|. Here we would like to focus on the different 
situation (Fig. [1]) considering the phonon modes coupled 
to the side-attached quantum dot. This is reminiscent of 
the setup discussed by M. Biittiker [11] except that the 
fermion reservoir is here replaced by the phonon bath. 

In analogy to (O we express the selfenergy ^^(ti^) of 



QD2 by the following contributions 

S°(w) =tGi(cj)r + sf (. 



(12) 



where tGi{uj)t* originates from the interdot hybridiza- 
tion and the second term is induced by the phonon reser- 
voir. Since QD2 is assumed to be weakly coupled with 
the central dot we approximate the selfenergy 1!i^{uj) 
adopting the local solution. The selfenergy H^(uj) can 
be determined by means of the Lan g-Fi rsov canonical 
transformation which effectively gives j29{ 



(13) 



£2 



Sf(a^) 




with the quasiparticle energies £2(0 = ^2 — Ap -I- ZllIq, the 
corresponding polaronic shift Ap = A^/wq and tempera- 
ture dependent spectral weights A{1) [25, 2^ 



A{1) = e^^ ^nuo/keT 



As usually, we introduce a dimensionless adiabadicity pa- 
rameter g = Nphil) = [e'^o/'^bT _ -g ^-j^^ gQg^_ 

Einstein distribution and // denote the modified Bessel 
functions. In particular, for the ground state the equa- 
tion simplifies to 



lim A(l) = e-s^ 9(1). 
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FIG. 2: (color online) Spectral function of the side-attached 
quantum dot obtained at low temperature for the model pa- 
rameters £2 = — IFjv, U2 — SFiv, cjo = 0.2rjv, g = 1, ei — 0, 
A = IOFat assuming a weak interdot coupling t — 0.2rjv. 



In figure [5] we illustrate the characteristic spectrum 
of the side-attached quantum dot (see section IV. A for 
technical details). We notice two groups of the narrow 
peaks. The lower one starts from the energy 62 ^ s — Ap 
(where Ap = /ujq is the polaronic shift) followed by a 
number of equidistant phonon peaks spaced by wq. The 
upper phonon branch is sep arated by C/2 and it mani- 
fests the charging effect '2&\. For the coupling g = 1 we 
observe only about five phonon peaks but in the antia- 
diabatic regime (5 3> 1) their number considerably in- 
creases. Such tendency is shown in section III discussing 
the spectrum of QDi. Let us also stress that the in- 
terference peaks have a rather tiny but yet finite width 



C. Subgap transport 

Charge transport in a subgap regime \eV\ < A is gener- 
ated only by the Andreev mechanism. Electrons coming 
from the metallic lead are then converted into the Cooper 
pairs in superconductor simultaneously reflecting holes 
back to the normal lead. Such anomalous current lAiV) 
can be expressed by the Landauer-type formula [30j 



Ia{V) 



2e 



dojTAiiu) [f{u~eV,T)^f{ij + eV,T)]il6) 



where f{uj, T) = l/ [e'^/'^^T ^ ^^^^ Fermi-Dhac func- 
tion. Andreev transmittance 7^(1^) depends on the off- 
diagonal part of the Green's function Gi (w) via [s^] 



TA{io) = Tl |Gi,i2(c^)r 



(17) 



(15) 



The transmittance ([TT]) can be regarded as a qualitative 
measure of the proximity induced on-dot pairing. Under 
optimal conditions it approache s unity when uj is close 
to the quasiparticle energies zb-^/e^ + (Ts/'^)'^- The An- 
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dreev transmittance is also sensitive to other struc- 

tures, for instance originating from the interferometric 
effects [Ullllll. 

In what foUwos we shall study the differential Andreev 
conductance 



Pdi(») 



Ga{V) 



dV 



(18) 



exploring its dependence on the phonon modes. We start 
this analysis assuming that both quantum dots uncorre- 
lated and we next extend it (in section IV) by considering 
the correlation effects. As a general remark, let us no- 
tice that the even function Ta(— w) = T4(a;) implies the 
symmetric Andreev conductance Ga{—V) = Ga{V) re- 
gardless of any particular features due to interference, 
phonons, correlations or whatever. Physically this is 
caused by the fact that particle and hole degrees of free- 
dom participate equally in the Andreev scattering. 



III. UNCORRELATED QUANTUM DOTS 

Upon neglecting the correlation selfenergies E^(cj) = 
one has to solve the following coupled equations 



G^\uj) = g^\..)-|i|^Gi(c.)-Ef (c.) 



(20) 



We have computed numerically the matrix Green's func- 
tions Gi{uj) for a mesh of energy points appropriate for 
the model parameters (mainly dependent on wo, 5 and t). 
Practically already about ten iterations proved to yield 
a fairly convergent solution. 

Since the proximity induced on-dot pairing predomi- 
nantly affects the energy region around we first con- 
sider the instructive case ei = 0, £2 7^ £i- In figure |3] we 
present the equilibrium spectrum pdi{uj) of the central 
quantum dot obtained for three representative coupling 
constants g = A/wq corresponding to the adiabatic limit 
g <^ 1 (upper panel), the antiadiabatic regime g ^ \ 
(bottom panel) and the intermediate case (middle panel) . 
On top of two Lorentzian peaks centered at the quasipar- 
ticle energies ±r5/2 we clearly see formation of the Fano 
resonances. They appear at energies e^ + Iojq and at their 
mirror refiections (because of the particle - hole mixing 
21|). Number of these phonon features depends on the 
adiabadicity parameter g. For the adiabatic regime there 
appear only a few phonon features whereas in the oppo- 
site antiadiabatic limit there is a whole bunch of such 
narrow structures. In the latter case they seem to have 
an irregular structure, but after closer inspection we can 
clearly see the Fano- type shapes (see the inset). 

Phonon driven replicas of the Fano lineshapes appear 
also in the Andreev transmittance (see Fig. |4|). In a dis- 
tinction to the spectral function pdi{uj) the resonances 
show up in a symmetrized way due to reasons mentioned 




FIG. 3: (color online) The interferometric Fano-type line- 
shapes appearing at T = in the spectral function pdiij-o) 
of QDi. Numerical calculations have been done for the un- 
correlated quantum dots Ui = Q using the model parameters 
£1 = 0, £2 = IFjv, t = Q.2Vn, Vs = SFat and A = IOFat. 



in the preceding section. Again we notice the broad max- 
ima centered at the subgap quasiparticle energies ±Fs/2 
accompanied by a number of the Fano-type resonances at 
± (£2 + ^i^o) ■ Spectroscopic measurements of the Andreev 
conductance would thus be able to detect such phonon 
induced interferometric features. 

The subgap quasiparticle lorentzian peaks (often re- 
ferred as the hound Andreev states) depend on the en- 
ergy level £1. In the case of single quantum dot (i.e. 
for vanishing t) the spectral function pdi{uj) consists of 
two lorentzians at ±i?i (where Ei — + (rs/2)^) 
broadened by F^r. Their spectral weights are given by 
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FIG. 4: (color online) The Andreev transmittance (upper 
panel) and the differential subgap conductance GAiV) (bot- 
tom panel) obtained for the intermediate coupling limit g = 1 
using the same parameters as in figure [3] 



the BCS factors ^ (1 ± ei/£'i). This fact has some im- 
portance also for the interferometric features. Figure [5] 
shows the spectrum pdi{^^) for several values of £i. When 
the energy £i moves away from the Fermi level (by ap- 
plying the gate voltage) we observe a gradual redistri- 
bution of the quasiparticle spectral weights accompanied 
with suppression of the Fano resonances, especially at 
— (£2 + 'wo). The Andreev transmittance Ta{oj) and dif- 
ferential conductance Ga (Y) are even functions therefore 
such particle-hole redistribution is not pronounced, nev- 
ertheless suppression of the phonon induced Fano line- 
shapes is well noticeable. 



IV. CORRELATION EFFECTS 



In this section we address qualitative effects caused by 
the Coulomb repulsion between the opposite spin elec- 
trons. Roughly speaking, we expect some possible sig- 
natures of the charging effect (Coulomb blockade) and 
eventual hallmarks of the Kondo physics. Since the elec- 
tron transport occurs in our setup via QDi we suspect 
that predominantly the Coulomb potential Ui can have a 
significant role. For completeness we shall however study 
the influence of correlations on both quantum dots. 




E, /r, 



FIG. 5: (color online) Dependence of the spectral function 
Pdi(aj) and the Andreev conductance Ga{V) on ei (tunable 
by the gate voltage). Calculations have been done for the 
same parameters as in figure [4] 



A. Influence of U2 

We start by considering the effects of finite U2, ne- 
glecting correlations on the central quantum dot (Ui =0). 
Since the side-attached quantum dot is weakly hybridized 
with QDi therefore the indirect influence of external 
leads on QD2 should be rather meaningless. For this 
reason we impose a diagonal structure of the selfenergy 



diag / 






-^diag 



(21) 



and here i = 2. We next approximate the diagonal terms 
of pT|) by the atomic limit solution 



1 



1 - "2,, 



(22) 



where | =t, f It has been pointed out [3l| that 

the selfenergy defined in equation (|22p coincides with the 
second order perturbation formula 



2.a (^) 



U2 n2.s~ 



{U2) "2,g(l"»2,g) 
W - £2 - C/2(l - n2.s) 



(23) 



and it can be generalized into more sophisticated treat- 
ments in the scheme of iterative perturbative theory [s^l . 
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FIG. 6: (color online) Spectral function pdi(ijj) of the central 
quantum dot (top panel) and the differential Andreev con- 
ductance (bottom panel) obtained for ei = 0, £2 = — IFat, 
t = 0.2riv, g ^ 1, uo ^ 0.2rjv, A = lOFiv assuming the 
Coulomb potential U2 = 5Tm- The corresponding Pd2{ijj) is 
shown in Fig. (2] 



For the weak interdot coupling t we expect however that 
corrections to (|22I23I) are not crucial. We skip here the 
higher order superexchange mechanism leadind to the ex- 
otic Kondo effect 33] which is beyond the scope of our 
present study. 

The top panel of figure [5] shows the spectrum pdi{uj) 
obtained at low temperature for t/2 = STjv, C^i = 0. As 
far as the side attached quantum dot spectrum is con- 
cerned it reveals a bunch of phonon peaks formed near 
the energy 62 and another group of states around the 
Coulomb satellite £2 + 1^2 (see figure [2]). These phonon 
signatures appear in pdii^^) as the Fano-type resonances. 
Due to the absorbed superconducting order we can no- 
tice effectively four groups of such Fano-type structures 
nearby the energies £2, £2+C^2 and at their mirror reflec- 
tions. The Andreev transmittance Ta{i-o) is symmetrized 
versions of what is shown in figure [S] therefore the result- 
ing differential conductance is even function of applied 
voltage V (see the bottom panel in Fig. [2]). 



B. Influence of Ui 

Correlations originating from the Coulomb repulsion 
Ui have a totally different effect on the transport prop- 
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FIG. 7: (color online) Spectral function pdi{oj) of the cen- 
tral quantum dot obtained in the Kondo regime for ksT = 
10-^Fiv, ei = -2F]v, U2 = ISFjv, Fs = 4Fiv, t = 0.2Fjv, 
Ldo = 0.2Fjv, g — 1, A ^ Tn. The vertical arrows indicate 
positions of the subgap Andreev states at ±£1 and at their 
Coulomb satellites. We can notice that the Kondo peak at 
cj = is distinct from the phonon induced Fano resonances. 



erties than above discussed U2- The central quantum dot 
is directly coupled to both external leads therefore on one 
hand by it experiences the induced on-dot pairing (due to 
Fs) and, the other hand, the Kondo effect (due to Fat). 
These phenomena are known to be antagonistic. Their 
nontrivial competition in a context of the quantum dots 
has been discussed theoreti cally by many groups using 
various techniques (see Ref. [1^ for a survey). 

To recover basic qualitative features we shall follow 
here our previous studies [s^ which proved to yield sat- 
isfactory results for the single quantum dot on interface 
between the metallic and superconducting leads [l^ . We 
choose the correlation selfenergy E^(a;) in the form (pij) 
and determine its diagonal parts by the equation of mo- 
tion approach [35| . Formally we use 

sf7(u;) = C/iK,^-Ei(u;)] (24) 
Ui [ni,g-S]i(a;)] [£3(6^) + t/i(l-ni,^)] 
w - £1 - So(w) - [S3(c^) + Ui{l- ni,ff)] ' 



where 



k 



1 



1 



/(Ckiv) 
1 



t^-CkW CJ-C/i- 2£i-|-^kA' 

for i> — I 
for ly = 3 



(25) 



and as usually i;o('^) = Sk l^k7Vp/(w-CkAr)==-»F7v/2. 

Let us remark that upon neglecting the terms Si(aj) 
and 23(0;) the selfenergy ((M)) nearly coincides with the 
second order perturbation formula (j23p 



lim 



Ediaq / \ 



(26) 
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FIG. 8: (color online) The differential Andreev conductance 
Ga{V) for tlie set of parameters discussed in figure [7] and 
several values of £2- 



except of I]o(w) present in the numerator of (|27p. This in- 
dicates that equation ([M)) is able to capture the charging 
effect (Coulomb blockade). The additional terms S,/(w) 
provide corrections which are important in the Kondo 
regime, i.e. for ei < < £1 + C/i at temperatures below 

■n-ei(ei + C/i) 



Under these con- 



Tk = 0.29^f7irA,/2 exp 

ditions at there forms the narrow Kondo resonance 
of a width scaled by ksTK- In our present method such 
Kondo resonance is only qualitatively reproduced (spl- 
its structure in the low energy regime < ksTx must 
be inferred from the renormalization group or other more 
sophisticated treatments. 

Let us now point out the main properties character- 
istic for the Kondo regime. The equilibrium spectrum 
of QDi illustrated in figure [7] consists of four Andreev 
bound states (indicated by the vertical arrows) centered 
at ±^ei + (rs/2)2 and ± ^(ei + Uif + [Vs/if. The 
fact that El is located aside the superconducting en- 
ergy gap causes asymmetry of the quasipartice spectral 
weights. Besides these broad lorentzians we additionally 
notice the narrow peak at the Fermi level signifying the 
Kondo effect. Such Kondo peak is considerably reduced 
in comparison to the normal case Fg = because of a 
competition with the on-dot pairing [3J| . On top of this 
picture we recognize the phonon degrees of freedom ap- 
pearing as the Fano-type resonances at ±(^2 + l^o)- 

The above listed effects are also detectable in the dif- 
ferential Andreev conductance (see Fig. [5]). Ga{V) has 
local maxima at voltages corresponding to the energies 
of the subgap bound states. Furthermore, similarly to 
our previous studies js^l , we notice that the Kondo peak 
leads enhances the zero-bias Andreev conductance. This 
property has been indeed observed experimentally p^ . 
In the present situation we additionally observe the Fano- 
type resonances. They destructively affect the zero-bias 
enhancement whenever the phonon features happen to be 



located nearby the Kondo peak. The zero-bias feature it- 
self is also quite sensitive to the asymmetry ratio F5 /F^v 
- practically it is visible only when both couplings T p are 
comparable [sij . 



V. CONCLUSIONS 

We have studied influence of the phonon modes on 
the spectral and on the Andreev transport in the double 
quantum dot vertically coupled between the metallic and 
superconducting electrodes. Our studies focused on the 
weak interdot coupling t assuming that phonons directly 
affect only the side-attached quantum dot. Under such 
circumstances an external phonon bath leads to some in- 
terferometric effects, reminiscent of the dephasing setup 
introduced by Biittiker psj . 

In particular we find a number of the equidistant Fano- 
type patterns manifested both in the effective spectrum 
and in the subgap transport properties. These lineshapes 
appear at ±(£2 + ''^o) (where wo is the phonon energy, I 
is an integer number and £2 = £2 ~ A^/cjq describes the 
QD2 energy shifted by a polaronic term) . They can be re- 
garded as replicas of the initial interferometric structures 
in absence of the phonon bath formed at ±£2 [20, 23]. 

Electron correlations Ui on the quantum dots can in- 
duce additional Coulomb satellites of these Fano features. 
We have investigated in some detail how the correla- 
tion effects get along with the Fano interference tak- 
ing into account the induced on-dot pairing. We notice 
that phonon features are sensitive to the subgap Andreev 
states (dependent on Ui) and to the Kondo effect. The 
latter one is important whenever the phonon lineshapes 
are induced in a vicinity of the Kondo peak. We thus 
expect that quantum interference would destructively af- 
fect the Kondo physics by partly suppressing its zero-bias 
hallmark [l^. 

It would be of interest for the future studies to check if 
the presently discussed effects are still preserved when 
the interdot coupling t is comparable to the external 
hybridization F^. We suspect, that the Fano-type pat- 
terns shall evolve into some new qualities typical for the 
complex molecular structures. Furthermore, the role of 
Coulomb interaction U2 might prove to be more influen- 
tial via the higher order exchange integrations inducing 
some exotic kinds of the indirect Kondo effect 33 1 . These 



nontrivial issues deserve further studies eventually using 
some complementary methods. 
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Appendix A: Fano-type interference: 
phenomenological arguments 

Here we would like to explain in a simple way why 
the interferometric Fano structures show up in the trans- 
port properties of DQD system. In general, the Fano- 
type lineshapes [11] emerge whenever the localized (res- 
onant) electron waves interfere with a continuum (or 
with sufficiently broad electron states). Following [33| 
let us consider the very instructive example in which 
the "direct" transmission channel td — VGd.e'^'^'^ (here 
Gd denotes its conductance and (pd stands for an ar- 
bitrary phase) is combined with the transmission am- 
plitude tr(uj) = JWr ^^i'tr^"ip^ ^ /o of another "res- 

onant" level Er- From general considerations [35| the 
corresponding conductance of such "resonant" level is 
Gr = ■ Effectively these two channels yield 

the following asymmetric structure 



1 ' 



(Al) 



Similar reasoning can be applied to the T-shape dou- 
ble quantum dot system shown in Fig. [1] For simplic- 
ity let as neglect the phonon bath and assume that 
both electrodes are normal conductors. In the case 
of weak interdot coupling <^ r| the side-attached 
dot QD2 plays the role of "resonant" channel with its 
transmission amplitude t^(cj) ~ v^r -^-tttttat , where 

Gr 



' w-e2-Hit/2 ' 

=j^. On the other hand the other central dot 
QDi provides a relatively broad background td{uj) = 

2e2 4rjvrs (riv+rs)/2 p onpro-lpQ / 1 ^ <r„ 

IT {T«+Ts? u^-e^+i(TM+Ts)/2- energies a; - £2 

the latter amplitude is nearly constant td{ijj) — \/Gde^'^^ 

1 



with y/Gd = 



2e2 



4rwrs 



2(ei-e2)/(rN+rs)+i 



Un- 



h (rjv+rs)^ 

der such conditions the resulting conductance G{uj) = 
\td + tr{^)\ indeed reduces to the Fano structure (lAip . 
Some more specific microscopic arguments in support 
for the Fano-type interference of the strongly correlated 
quantum dots have been discussed at length e.g. by 
Maruyama [13] and by Zitko [ll]. 



where tj = (cj — e,.)/(ir) and q = i 
characteristic asymmetry factor. 
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